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Abstract 

A particle produced in a hard collision can lose energy through bremsstrahlung. It has long been 
of interest to calculate the effect on bremsstrahlung if the particle is produced inside a finite-size 
QCD medium such as a quark-gluon plasma. For the case of very high-energy particles travel- 
ing through the background of a weakly-coupled quark-gluon plasma, it is known how to reduce 
this problem to an equivalent problem in non-relativistic two-dimensional quantum mechanics. 
Analytic solutions, however, have always resorted to further approximations. One is a harmonic 
oscillator approximation to the corresponding quantum mechanics problem, which is appropriate 
for sufficiently thick media. Another is to formally treat the particle as having only a single sig- 
nificant scattering from the plasma (known as the A^=l term of the opacity expansion), which is 
appropriate for sufficiently thin media. In a broad range of intermediate cases, these two very dif- 
ferent approximations give surprisingly similar but slightly differing results if one works to leading 
logarithmic order in the particle energy, and there has been confusion about the range of validity 
of each approximation. In this paper, I sort out in detail the parametric range of validity of these 
two approximations at leading logarithmic order. For simplicity, I study the problem for small 
and large logarithms but aglog <C 1. 



I. INTRODUCTION AND RESULTS 



A. Background 

There is a prototypical toy problem often considered in theoretical discussions of gluon 

bremsstrahlung in a QCD medium such as a quark-gluon plasma: Consider a high-energy 
quark or gluon that is produced by some hard scattering event and then propagates through 
a length L of a uniform QCD medium before emerging into vacuum. What is the effect of the 
medium on the probability for gluon bremsstrahlung from this high-energy particle? This 
is known as the brick problem.^ The problem is complicated by the LandaTi-Pomeranchuk- 
Migdal (LPM) effect [1, 2]. The quantum mechanical duration (formation time) of the 
bremsstrahlung process grows with increasing energy and eventually exceeds the mean free 
time between collisions. As a result, successive collisions of the high-energy particle with 
the plasma cannot be treated as independent from each other for the purpose of calculating 
the probability of bremsstrahlung. 

There is a general formalism for treating this problem,^ but analytic solutions have re- 
quired additional approximations. Baier, Dokshitzer, Mueller, and Schiff (BDMS)'^ [14] in- 
vestigated the problem in the limit that the energy was high enough, and the medium thick 
enough, that the number of collisions within the bremsstrahlung formation time was large 
— so large that InA^ could be treated as large. In this limit, they made an approximation, 
known as the harmonic oscillator (HO) approximation, that reduced the general formalism 
to a certain type of harmonic oscillator problem. They solved for the medium effects on 
the spectrum of gluon bremsstrahlung, to be reviewed below. Prom the spectrum, they 
computed the size AE of the medium effect on the average energy loss of a high-energy 
particle of energy E.^ The qualitative form of their result depends on the thickness L of 
the medium compared to the typical formation length for gluon bremsstrahlung in an 



^ There are other versions of this problem. Sometimes people consider the case of a high-energy quark or 
gluon that propagates a relatively long distance through vacuum, then enters a uniform QCD medium of 
length L, passes through it, and exits the other side, approximately maintaining its direction throughout. 
I instead consider the case where the particle is first created inside the medium. Created could mean as 
one member of a particle/anti-particle pair well separated in angle (in which case one would separately 
compute the medium effect on bremsstrahlung from the other particle), or it could mean the final state of 
a large-angle deflection of a pre-existing particle (in which case one would separately compute the medium 
effect on initial-state radiation). 

^ See [3-8] for the original development. See also [9] for a summary in a language that generalizes naturally to 
the problem of non-fixed scatterers, and for a discussion of how the formalism is related to that developed 
in Refs. [10-12] for the case of infinite media. For a nearly complete calculation of bremsstrahlung in the 
infinite mcdinm case, to leading order in a, see Ref. [13]. 
See also the earlier work with Peigne of Refs. [5, 6]. 

* The single number given by the average energy loss leaves much to be desired as a description of the 
final-energy probability distribution because that distribution tends to have large, non-Gaussian tails. 
See the discussion in Sec. 3 of Ref. [15] or Ref. [13]. However, here my purpose is just to use it as an 
example for the sake of theoretically comparing the roles of the HO and A'^=l approximations. 
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infinite medium, which is parametrically 

i^oo~y^. (1.1) 

Here, q is the typical squared transverse momentum per unit length transferred via elastic 
collisions to a high-energy particle as it traverses the medium (more discussion later). For 
thick media (L ^ L^o), they found that AE grows linearly with L, as one would expect. 
For thin media (L -C -^oo), they found [14]^ 

AEho ^ ICMaL' ^ TrCC-AaWL^ In (^) , (1.2) 

to leading order in inverse powers of the logarithm. Here s is the species (quark or gluon) 
of the high-energy particle, and Cr is the quadratic Casimir of a given color representation. 
A/" is the density n of plasma particles weighted by group factors as® 

^f=j^ [Ca% + CF(nq + riq)] = I rig + |(nq + n^) = (l + iiVf)T^ (1.3) 

where is the number of quark flavors. 

In contrast, various other authors have investigated the opposite approximation, starting 
from early work by Wiedemann and Gyulassy [16] and by Gyulassy, Levai, and Vitev (GLV) 
[17-19]. Instead of treating the number of elastic collisions as large, they expand order 
by order in the number of collisions. This is known as the opacity expansion. The leading 
term, corresponding to iV = 1, gives^ 

AEn=i ~ ttC^Ca^WL^ In ( ^-r] (1.4) 



^ Specifically, the first equality in (1.2) is equivalent to Eq. (49) of Ref. [14], which can be expressed in 
terms of q as {AE)/L = jaCx^sL, where Ca = Nc for SU(A^c) gauge theory. Since qs is proportional 

to Cg, one can rewrite this in the form jaCgqAL, which will be more convenient for my later discussion. 
The last equality in (1.2) is given by the formula for which I review later in (1.13b). 
^ Here rfn is the dimension of color representation i?, rig is the total gluon density, and riq and riq are the 
total quark and anti-quark densities, summed over flavor. One may equivalently write J\f = (^A/dA)?^g -l- 
{tF/dp){Ti(^ + fiq) — 2tp^n+ + ANftpTi- where Ir is the trace normalization defined in terms of color 
generators by tr(T^T]^) = t_R(5"'', and n± = J {2n)~^d^p {e^^ =F 1)~^ is the number density of a single, 
massless, bosonic/fermionic degree of freedom. 

Eq. (16) of Ref. [18] gives {CrolL'^h'^ /A\^\n{E/ii), where Ag is the gluon mean free path and /x"^ is 
the color electric screening length. The older literature on gluon bremsstrahlung often unnecessarily 

normalizes answers in terms of Ag, which is not well defined in leading-order perturbation theory because 
of a logarithmic infrared divergence from magnetic scattering. But the result for bremsstrahlung does 
not depend on these details. Using the purely electric scattering models assumed in older calculations, 
/i^/Ag = AnCp^a^M. This substitution recovers independence from the details of electric vs. magnetic 
screening. The fact that the appropriate lower scale in the logarithm is of order miL can be found in 
the work of Zakharov [20] and is nicely laid out in the presentation of Salgado and Wiedemann [21]. 
Alternatively, readers of GLV can see it by noting that Eq. (15) of Ref. [18] (Eq. (130) of Ref. [19]) is only 
valid when 7 <C 1 and so x ^ Lml/E, and so the infrared logarithmic divergence in the x integration of 
that equation is cut off by this lower bound on x and generates the logarithm in (1.4) above. 
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in the high-energy hmit, to leading order in inverse powers of InlE/rn^L). 

The results (1.2) and (1.4) from opposite assumptions about the relevant number of col- 
lisions are surprisingly similar, differing only in the argument of the logarithm. Two natural 
questions arise. Which formula is correct for what range of media thickness LI Which 
description captures the correct physics: Is there a single collision with the medium which 
dominates the medium's contribution to bremsstrahlung energy loss, or is AE dominated 
by processes where many scatterings are important? 

The quahtative difference between the HO and A^=l approximations becomes more pro- 
nounced if one looks more generally at the gluon bremsstrahlung spectrum instead of focusing 
on the single number AE. For the brick problem, the HO approximation gives the result 
[14]8 

d (X 
OJ - i"vac)HO = -xPs-^g{x) In |cos(a;oL)| (1.5) 

with 

2 . [(1 - x)qA + x'^Qs] 

LOr, — —I — . (I.D) 

° 2x{l-x)E ^ ' 

Here, / is the probability of gluon bremsstrahlung, with the corresponding probabil- 
ity if the process which created the high-energy particle had instead taken place in vac- 
uum. Ps^g{x) is the usual vacuum splitting function,^ uj is the energy of the (high-energy) 
bremsstrahlung gluon, and x = uo/E is lis momentum fraction. |a;o|~^ is of order the for- 
mation length /oo(^) for a bremsstrahlung gluon of frequency uj in an infinite medium. The 
previous formula (1.2) for AE is just the uj integral of (1.5) in the limit L -C L^o- In that 
limit, the uj integral is dominated by small x such that loo{oj) ~ L. But now fix uj and 
consider thinner and thinner media such that L is small compared to the formation length 
loo{oj). The small L limit of (1.5) is 

d ,r T \ n / \ \oJoLy^ ,^ 

OJ —{I - /vac)HO ^ -xPs^si^) (1.7) 
acj TT iZ 

Focusing on the case x <^ 1, ior simplicity, gives 

c.^(/-/vac)HO-2^L — Lin ^-^j, (1.8) 

where the last equality relies on the formula for q, to be reviewed momentarily. In contrast, 
the N=l term of the opacity expansion gives 

U!—{I - Iy^)N=l^ L. (1.9) 

duj UJ 

Because it is proportional to rather than L"^ . the A^=l result clearly dominates over the 
HO result for small L. The HO result seems completely at odds with the A^=l result for 
small L. When is the HO result correct? 

All of these issues were raised some years ago by Zakharov [20] . He concluded that the 
HO analysis should not be trusted in cases when the medium thickness L is less than or 



For a relatively simple formula for more general situations of expanding, inhomogeneous media, see Ref. 
[9]. See also the earlier work of Ref. [22]. My sign convention in (1.6) is that of Ref. [9]. 
Pq^g(a;) = Cf[1 + (1 - xf\lx; Pg^g(a;) = Ca[1 + a;^ + (1 - xf\lx{\ - x). 
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order the relevant formation length Loo or Zoo(<^)- In this paper, I return to this problem and 
show that the HO approximation is valid over a wider range of L, and I elucidate in more 
detail the interplay between contributions to bremsstrahlung (i) arising from large numbers 
of scatterings and (ii) dominated by a single scattering. 

I will assume that the particle energy E is so large that \n{E/T) can be treated as a 
large number, where T is the plasma temperature. Because I want to pursue a qualitative 
understanding of how the HO and A^=l results fit together, I will usually just focus on the 
parametric form of formulas. Though I will treat logarithms as large, I will formally assume 
that a is so small that aln{E/T) is small. So, for instance, I will ignore running of the 
coupling and treat a as fixed. The purpose of these various limits is to provide a clean, 
theoretical situation for conceptually disentangling the HO and A^=l approximations. For 
experimentally achievable quark-gluon plasmas, of course, logarithms are not huge and a is 
not tiny.^^ 



B. Transverse momentum diffusion 

To describe my results, it is useful to first characterize the total transverse momentum 
Qj. that a high-energy particle picks up, due to screened Coulomb-likc interactions, as it 
crosses length L of a QCD medium. In a perturbative quark-gluon plasma, the differential 
elastic scattering rate for a high-energy particle to pick up transverse momentum q± is 

for q± ^ T. The behavior is similar, with a slightly different coefficient, for smaller q± 
down to mo, where Debye screening kicks in. Over multiple collisions, the net transverse 
momentum transfer will random walk, and its average will be 

(Qi)avg d'q^f^ql (1.11) 

Using (1.10), this integral has a well-known logarithmic divergence which is cut off in the 
infrared by Debye screening. The UV end is cut off by the kinematic limit g™^^ ~ V ET, 
where the plasma temperature T gives the typical energy of a plasma particle. However, the 
average Q\ will not be an interesting quantity for our purposes. The probability that there 
is at least one collision with individual momentum transfer of order q± over the distance L 
is of order 



So it is unlikely to have any individual collisions with q\ 3> Cro?ML. If we are interested 
in the typical (i.e. median) Q\ instead of the average Q\, we should use this value of as 



For some discussion of running coupling in the bremsstrahlung problem, see for example Sec. VI of Ref. 
[24], which combines earlier observations of Refs. [6] and [23]. 

There is some theoretical information concerning the efficacy of expansions in 1/ lii{E/T) in the context 
of infinite-medium bremsstrahlung calculations for weakly-coupled plasmas. In Ref. [24], it was found 
that the error of making a next-to-leading logarithm approximation is < 20% when E > 10 T. 
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FIG. 1: Probability distribution for a high-energy particle to pick up total transverse momentum 
via 2 — > 2 scattering while traveling through a medium of thickness ~ L. All formulas show 
only parametric dependence. The dashed line shows the behavior of the HO approximation at 
large Q±, and the dotted line shows the behavior of the single-scattering {N=l) approximation at 
low Q±. The kinematic upper limit Q\ ~ ET of the high-Qx tail is not shown. 



an upper cut-off on the integration in (1.11).^^ The result is then 

(Qi)typ = QrL (1.13a) 

with^^ 

QR ^ AnCna'Afln (^^^^) (l-13b) 

for a particle with color representation R. Throughout this paper, I will use q to denote 
the typical (rather than average) squared transverse momentum acquired per unit length, as 
given by (1.13b). 

A simple way that one can get to an equivalent result is to self-consistently use Q± itself 
to cut off the UV logarithmic divergence in (1.11), so that 

(gi)typ AnCRa'N-Lln (iSl^) . (1.14) 

This gives the same result as (1.13) up to corrections that are subleading in inverse powers 
of the logarithm. We can also write it in the form 

(gi)typ ^ ^nCRa'AfLln (^) . (1.15) 



This important distinction of the typical or characteristic Q'^ as opposed to the average was made previ- 
ously in Sec. 3.1 of Ref. [6]. 

In the context of infinite media, see Refs. [25, 26] for the weak-coupling evaluation of q beyond leading log 

order and its application to bremsstrahlung calculations. See also the related, earlier work of Ref. [27]. 
Throughout this paper, I will treat large logarithms as parametrically large, but I will treat logarithms of 
logarithms, such as ln\n.{Q\/m^), as being of order 1. 
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A qualitative sketch of the probabihty distribution of total is shown by the solid line 
in Fig. 1. Typical events show a Gaussian peak characteristic of diffusion in Q± space, 
whose width is given by (1.13). However, there is also a large Q± tail of rare events, where 
one of the collisions with the medium has q± ^ {Q±)typ- The probability distribution for 
these events is simply given by dTf,\/dq\ times L. The formulas shown in the figure are all 
parametric and do not show multiplicative factors of 0(1). They also do not show group 
factors such as Cr. The normalization 1/qL of the height of the diffusion peak can be 
determined from the requirement that the total probability is 1. By comparing the single 
scattering formula of the high-(5± tail and the Gaussian formula for a diffusion peak, one 
can parametrically estimate that the transition between the two occurs when the probability 
distribution is down from the peak value by a factor of order 

I will use the short-hand notation "log" defined above to denote this particular logarithm in 
figures. Further review of the important aspects of Fig. 1 is given in appendix A for readers 
desiring more detailed explanation. 

The HO approximation corresponds to ignoring the large-Qx tail of this distribution and 
approximating the probability distribution as a standard diffusion Gaussian peak. 



scatt 



HO lV±;typ 



'typ. 



(1.17) 



depicted qualitatively by the dashed line in the figure. The A'"=l approximation, in contrast, 
involves using the single scattering formula oc l/Q^ ^r all momenta, all the way down to 
the Debye mass. This is depicted by the dotted line in the figure, but the low- momentum 
cut-off at Q\ ~ ml is not shown. The double arrows beneath the plot indicate over which 
regions these two approximations are good approximations to the actual distribution. A 
cartoon of a typical scattering is shown in Fig. 2a. In contrast, a corresponding cartoon of 
one of the rare high-Qj. events is shown in Fig. 2b. Note that there are still many scatterings 
in this case, but a single one of those scatterings dominates Q±. I will assume throughout 
this paper that the medium is thick enough that the high-energy particle undergoes many 
soft collisions on its way through, corresponding to Fig. 2a or Fig. 2b. Parametrically, this 
assumption is that L ^ {C RaT)~^ . 

The probability of having more than one scattering with Qj^ ^ (Q±)typ is a parametrically 
small correction to Fig. 1 and so need not be considered. 

Throughout this paper, I will assume that energies are high enough that scattering and 
bremsstrahlung can be treated as nearly coUinear. In particular, I will restrict consideration 
to the case u) Q^. 



C. Results 

In this paper, I will show that the leading log result for /S.E gets two different types of 
contributions when L <^ L^: 



:i.i8) 
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FIG. 2: A depiction of (a) typical multiple scattering (the HO approximation) vs. (b) the rarer case 
(related to the A^=l approximation) where the total deflection is dominated by a single scattering 
with unusually large momentum transfer. The scattering angles are all exaggerated in this figure 
for the sake of visibility. The net angular deflection in both cases is parametrically ^ 1, and the 
small-angle scatterings are meant to be significantly smaller than the single larger (but still small) 
angle scattering in (b). 



The first logarithm is just the HO approximation of (1.2), corresponding to bremsstrahlung 
involving typical scattering from the medium. The second logarithm is due to events in- 
volving the rarer scatterings corresponding to the large tail of Fig. 1. Amusingly, the 
sum of these two logarithms simply gives the same mathematical formula as the full A^=l 
result of (1.4). However, depending on L, (1.18) can be dominated by the HO contribution. 
The formula (1.18) is qualitatively similar to a result by Zakharov [20],^^ but Zakharov's 
conclusion about the domain of applicability of the HO approximation was slightly different. 

If one thinks of taking the high-energy limit £" — > oo with fixed L, then the N=l term in 
(1.18) obviously dominates. But now consider fixing a large value of E and varying L. It's 
useful to rewrite (1-18) parametrically in terms of the typical formation length of (1.1): 



AE ~ CsCAa^NL^ 



villi V, 

V Loo J 



:i.i9) 



For L kept equal to any fixed fraction of Lqo, now the HO contribution dominates as — > oo. 
The two logarithms are equal when L is of order 

1/6 

L^. (1.20) 



rri 



D \ 

^aE) 

Note that is small compared to the typical formation length in the high energy limit. 
For L extremely large compared to L*, the result for AE will be dominated by the HO 



Specifically, see Eq. (23) of Ref. [20]. This result is only qualitatively similar to my (1.18) because of 
the upper limit Wcr on Zakharov's HO term, which in my notation he takes of order uj^r ~ qLi^ ■ The 
parametric treatment of this cut-off means that his HO contribution is only parametrically of order the 
HO result at leading-log order (which is all he asserts in his text). In my (1.18), the first term is exactly 
the HO result at leading- logarithm order. 
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approximation, corresponding to scatterings like Fig. 2a. In this limit, the formulas (1.2) 
and (1.4) in fact generate equivalent answers (that is, their difference is small compared 
to the result). For L extremely small compared to L*, will be will be dominated by 
bremsstrahlung involving the large-Qx tail of Fig. 1 and so by scattering like Fig. 2b. As 
I shall discuss, in this limit the physics of bremsstrahlung is effectively single scattering 
physics, even though there are multiple additional soft scatterings depicted in Fig. 2b. 

The parametric results for A£' are depicted quahtatively in Fig. 3, where the dotted hne 
denotes the full A^=l formula (1.4), and the double arrows below the graph again indicate 
whether the physics of the underlying elastic scattering is dominated by the HO or N=l 
type events of Fig. 2. One consequence is that the HO approximation remains valid when 
L ~ Loo- 

I should clarify that the scales of the axis in my figures are elastic and should not be 
interpreted as linear, though they do start at zero in the bottom-left corner. 

I will also preview my results concerning whether the bremsstrahlung gluon spectrum 
decreases as (the HO prediction) or (the N = 1 prediction) for small L and fixed 
gluon frequency uj. Consider 

APbrem = 0;^(/-/vac), (1-21) 

which parametrically is the medium effect on the probability of bremsstrahlung production 
of a gluon with frequency of order u). Fig. 4 gives a qualitative sketch of my result for APbrem 
versus the medium length L. In this figure, l^o is short-hand for the formation length Zoo(<^) 
for gluons of that frequency. For simplicity, I restrict attention to the case where 1—x is not 
small, in which case 

iM^xl^. (1.22) 

V Q'A 

Follow the L axis from right to left. As L drops below /oo, the curve follows the L"^ behavior 
of the HO approximation (1.8). At L ~ l^/^/\og, this corresponds to a drop of 1/log^ 
in probabihty from L ~ l^o- As L drops below Zoo/v^og, the behavior of the A'"=l 
approximation takes over. 

It's useful to note that one could generate the entire L behavior of Fig. 4 from the 
HO and A^=l results (1.8) and (1.9) if one made the assumption that the larger result is 
the correct one. 
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L 



'oo/\/log loo ioo/a 

JV=1 HO 



FIG. 4: Medium modification to probability for emitting a high-energy bremsstrahlung gluon of 
frequency ~ u: total probability vs. medium length L. 



Once again, one consequence is that the HO approximation remains valid when L is of 
order the relevant formation time, treating logarithms as large. However, in this particular 
case (unlike Ai?), the length scale at which the A^=l result takes over is smaller by only a 
square root of a logarithm. 

In the remainder of this paper, I derive and explain these results. In the next section, I 
will start with bremsstrahlung in QED rather than QCD plasmas. In particular, in section 
II A, I outline how the typical and rare scattering events of Figs. 2a and b can both give 
potentially important contributions to the total bremsstrahlung rate. It's then a matter of 
discovering which is the most important. Next I briefly review the scales associated with the 
LPM effect and then follow with detailed parametric estimates of the relative importance 
of the different cases. I move on to QCD in section HI, which requires relatively minor 
modifications to the QED analysis, although the final results for energy loss as a function 
of medium thickness are qualitatively quite different. The conceptually most important 
result in this development will be Fig. 20 (or Fig. 15 in the case of QED), which shows 
the relative importance of typical and rare scattering events to bremsstrahlung. Finally, in 
section IV, I reconcile the results of this paper with earlier analysis by Zakharov [20] . Some 
more detailed arguments concerning some of the qualitative points in this paper for QED 
and QCD bremsstrahlung are left to Appendix B and C, respectively. 

Throughout the main text, I will focus on gluon bremsstrahlung s — > gs with x < 1/2. 
For the case of g — > gg, that's everything because of the identity of the final state particles. 
For q gq, however, there is an additional contribution to energy loss from 1/2 < x < 1. 
The significance of this contribution depends on whether one considers AE to be the energy 
lost by the quark or the energy lost by the leading parton. The thin-media formulas (1.2) 
and (1.4) use the latter definition. In Appendix D, I spell out the details and explain the 
simple way in which Fig. 4 changes in the limit a; — > 1. 

Some of the parametric formulas I will derive in this paper were derived earlier by BDMS 
[15]. However, they did not keep careful track of logarithms when comparing HO and 
A'"=l contributions, which is important for a discussion of disentangling which is the most 
important at the order of leading logarithms. The situation is also one of many topics 
discussed in a recent mini-review by Peigne and Smilga [28], who discuss energy loss in both 
QED and QCD plasmas but do not keep track of logarithms. 
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(a) 



(b) 



FIG. 5: A non-relativistic example of the probability for soft bremsstrahlung radiation being unable 
to resolve details of charge particle tracks that are smaller than the photon wavelength. 



FIG. 6: The same processes as Fig. 5, now viewed from a highly boosted reference frame. The 
deflection angles are very small and the bremsstrahlung nearly collinear, but I have exaggerated 
the angles for the purpose of drawing the picture. 

II. BREMSSTRAHLUNG IN QED 

In this section, I focus on QED plasmas, deferring the treatment of QCD until section 



A. Basic Picture 

It's useful to first think about the case of soft bremsstrahlung, for which the charged 
particle can be approximated as classical. Imagine that Fig. 2 corresponds to possible 
particle tracks, and we want to estimate the bremsstrahlung probability. As a reminder of 
the origin of the LPM effect, first consider the case of a non-relativistic particle and recall 
that light cannot resolve features smaller than its wavelength. Thus, bremsstrahlung from 
the track in Fig. 5a will look the same as that from Fig. 5b, provided the wavelength is large 
compared to the distance scale over which the two trajectories behave differently. If one 
now Lorentz boosts this situation to extremely high energy, then the size of the region the 
photon cannot resolve will grow by a Lorentz factor and is now called the photon formation 
length, while the photon wavelength in that direction shrinks by a Lorentz factor. The 
photon therefore cannot resolve the difference between the situations of Fig. 6a and 6b. 

In the high energy case, photon bremsstrahlung will be nearly collinear with the charged 
particle, which can be understood as a result of the boost. In the ultra-relativistic limit, 
closer collinearity means larger formation lengths. Another useful mnemonic to keep in mind 
is that a photon emitted at angle 9 from an ultra-relativistic particle is not very sensitive to 
particle deflections small compared to 9. So the photon emission angle will be less than or 
order the net deflection angle of the charged particle within the formation length. 

Fig. 7 is somewhat similar to Fig. 6 but shows the case of propagation through a medium 
whose size is small compared to the formation length. The start of the particle trajectory 





III. 
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FIG. 7: The equivalence of bremsstrahlung (a) with and (b) without medium interactions for the 
case where the length L of medium traversed is small compared to the formation length. 

corresponds to whatever hard process (not shown) originally launched the high-energy par- 
ticle in its approximate direction of motion. Bremsstrahlung photons cannot resolve the 
difference between Fig. 7a and Fig. 7b, and so the medium does not have a significant effect 
on the probability of bremsstrahlung. 

Now consider the case of rarer collisions that involve a single larger-than-typical elastic 
scattering, as in Fig. 8a. This scattering can then affect bremsstrahlung radiation at larger 
than usual angles, corresponding to photons with smaller formation times. If the angle is 
large enough, so that the formation length in that particular case becomes smaller than the 
0{L) distance between the rare collision and the start of the particle's trajectory, then the 
photon can resolve the difference between Fig. 8a and 8b. The rare collision is therefore 
a second chance for bremsstrahlung, independent of the original event that produced the 
high-energy particle and without any LPM suppression. As far as this photon is concerned, 
the process is similar to the A^=l process shown in Fig. 8c. 

Because of the relatively small formation length, some readers may wonder if the addi- 
tional typical-ajigle scattering events in Fig. 8a can provide additional, distinct opportunities 
for bremsstrahlung, as depicted in Fig. 9, and so ruin the equivalence of Figs. 8a and 8c. This 
does not happen because the short-formation-length photons we have considered in the rare 
scattering case of Fig. 8 are emitted at angles large compared to the typical net scattering 
angle during one such formation length. The scatterings shown in the two additional ovals 
of Fig. 9, for example, do not produce significant bremsstrahlung radiation at such large 
angles. 

The upshot of this discussion is that scattering with larger than usual angles is rarer 
but, when it does happen, the probability for an associated bremsstrahlung photon is higher 
because there is less LPM suppression. Which type of process dominates the medium effect 
on bremsstrahlung depends on which of these opposing effects on probability is the most 
important. 

B. Review of LPM Effect 

For QED, one of the usual approaches to qualitative estimates of the formation length 
is the following:^^ two space-time points Xi and X2 on the charged particle's trajectory lie 



For a nice, very brief review, see, for example, the introduction of Ref. [29]. 
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FIG. 8: Like Fig. 7, but now (a) there is a rare scattering with larger than typical deflection 
angle, which affects bremsstrahlung at larger angles, which corresponds to smaller formation length 
(shaded ovals). Case (b) is no longer approximately (a), but case (c) is. 



same larger angle 



FIG. 9: As Fig. 8 but considering the possibility of similar photon production from other, typical 
collisions along the trajectory. This possibility is suppressed. 



within one formation length if the relative phase K^{Xi — X2)^ for photon emission from 
those two points is <C 1. If the particle is moving nearly linearly at close to the speed of 
light, this condition becomes |a;i— a;2|(l ~ cos^) -C 1, where 6 is the angle between the 
photon and the charged particle. Changing -C to ~ then qualitatively defines the formation 
length Zf, which for small 6 gives ulfO"^ ~ 1 and so 

~ (2.1) 

A more general way to the same result is to consider how off-shell in energy the intermediate 
particle line is in a simple bremsstrahlung diagram like Fig. 10, which is 

where P = p + k is the original momentum and the m are the effective finite-temperature 
masses of the particles. The formation time is the quantum mechanical duration of the 
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FIG. 10: One of the diagrams contributing to simple bremsstrahlung from a single scattering from 
the plasma. (The other important diagram is the one where the bremsstrahlung comes from the 
initial high-energy particle line.) 



off-shell state, /f ~ {SE) ^. If one ignores the masses, (2.2) can be rewritten in the form 

^ 6E pk[P^-!^Y x{l-x)Ee^' 

For X not too close to 1 (i.e. 1—x not small), this is the same parametric estimate as (2.1). 
As stated earlier, I will focus on x < 1/2 and ignore the case of small 1—x in the main text. 

Now, to set some scales, consider QED bremsstrahlung in an infinite medium, which is 
dominated by typical scattering events. The dominant photons are those whose angle 9 
relative to the charged particle is of order the net deflection angle AO of the charged particle 
from the scatterings it experiences during one formation time: Bremsstrahlung at larger 
angles 6 ^ A9 is suppressed. On the other hand, because of multiple collisions, the average 
angle 6 that the photon makes with the charged particle during a formation time cannot be 
smaller than order A6. So (2.1) becomes 

In an infinite medium (or any medium larger than the formation length), the typical deflec- 
tion angle A6 of the charged particle in one formation length loo is 

^AQ\ Q±oo {qJ-co) ^ /„ p-N 

(A^)oo ^ ^ — , (2.5) 

where Q±oo ~ (^^oo)^^^ is the transverse momentum the charged particle picks up over that 
distance. Combining (2.4) and (2.5), 

[e^ 

/ooM~W— . (2.6) 

V 'i^ 

Energy loss in an infinite medium is dominated by the case uj ^ E where the photon carries 
away a significant fraction of the particle's energy. In this case, the formation length becomes 
Lao ~ \/E/q, just like the QCD result quoted in (1.1). 

For the sake of easy reference, and for comparing and contrasting QED and QCD 
bremsstrahlung, I have collected in Table I some of the formulas described here and in 
section III. 
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QED QCD 



li 


1 1 


1 1 


formation length (general) 


ZjW ~ (^(Ae)2 




e 






characteristic bremsstrahlung angle 


(Qi)typ for L < Zoo 


qL 




typical momentum transfer 




IM 




infinite-medium formation length 


Loo ~ loo{E/2) 


V y 


V 9 


dominant l^o for energy loss 


Q 




gA^CAaWln(^) 


relevant Q\ per length 



TABLE I: Summary of various parametric formulas for the LPM effect in QED and QCD. These 
formulas assume 1—x is not small, and the QCD formulas ignore the difference between gA and qp 
in the case uj ^ E. 



C. Bremsstrahlung for a given 

Let be the total transverse momentum transferred to the charged particle as it tra- 
verses a medium of length L, and consider the case L <^ /oo(^) of a medium that is thin 
compared to the typical infinite- volume formation length (2.6) for photons with some fre- 
quency uj. In this section, I will focus on how the medium effect on the bremsstrahlung 
probability depends on Q^. 

1. LPM suppression in thin media 

Consider a charged, nearly massless particle that undergoes a single scattering in vacuum. 
The probability that this scattering will produce a photon with frequency of order uj (for 
any uj < E) is of order a times a coUinear logarithm. For the vacuum case, integration over 
bremsstrahlung frequencies gives rise to an additional, infrared logarithm in the probability. 
In considering medium effects on bremsstrahlung, I am going to delay integration over 
frequency until the end and will for now just consider the probability of emission of photons 
whose frequencies are of order some scale uj. The coUinear logarithm comes from integrating 
over photon directions that are very close to either the incoming or outgoing particle track 
in an isolated collision. It plays only a limited role in medium effects on bremsstrahlung. For 
the moment, I will ignore the coUinear logarithm and simply take a to be the additional cost 
in probability of emitting a photon of frequency ~ uj when there is a collision. The angle 
this photon makes with the incoming and outgoing particle tracks is of order the deflection 
angle of the particle track, with discussion of the possibility of more-nearly coUinear photons 
deferred to later discussion of the coUinear logarithm. 

In a medium, divide the multiple collisions along a trajectory into sets which are each 
roughly one formation length long. The photon cannot resolve the difference of a single 
vs. multiple collision in each set, but it does see each set as a distinct opportunity for 
bremsstrahlung. There is then an 0{a) probability for a photon emission from each set for 
which ~ 9. So, for instance, the shaded oval in Fig. 8a is associated with a probability of 
0{a) [times a coUinear logarithm] for emitting a bremsstrahlung photon at the angle shown. 

In Fig. 7a, there is also an 0{a) probabihty of emitting such a photon at the angle 
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Ql 



FIG. 11: Suppression factor e of medium effects as a function of Q\ for the case wliere the medium 
thickness L is smaller than the typical formation length looi}^) of the infinite-medium case. 



~ shown there. But the difference in emission probability with the vacuum case of 
Fig. 7b is small. In this case, I will write the medium contribution to the 0{a) probability 

of photon emission as 0{ea), where e is an LPM suppression factor due to the photon's 
failure to resolve (i) the collisions in the medium from (ii) the event that originally created 
the charged particle: 

A(cost of bremsstrahlung emission) ~ ea. (2.7) 

Fig. 11 shows the behavior of e as a function of Q\. I will explain this figure one feature at 
a time. 

For a given scattering trajectory — e.g. the typical scattering events of Fig. 7a or the 
rare events of 8a, or something in between — the relevant formation length will depend on 
the deflection angle, which for thin media will be related to the net transverse momentum 
transfer (5± while traversing the medium by 

~ %. (2.8) 

The corresponding formation length (2.1) is 

1 E'^ 

1 _. (2.9) 



Note that in a rare scattering case hke Fig. 8a, I should have estimated the formation 
length based on the net deflection over the formation length, shown by the shaded oval, 
rather then from the deflection over the entire trajectory. However, in this case, the rare 
scattering dominated the total angular deflection in any case, and so I do not need to 
distinguish between the two. 

The region e ~ 1 of Fig. 11 corresponds to cases k <^ L where the scattering process 
looks hke Fig. 8a rather than Fig. 7a. From (2.9), this condition is equivalent to 

Ql > — for e ~ 1. (2.10) 
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FIG. 12: Example of two scatterings with zero total deflection angle. 



2. The size of e 

Now turn to the case Q\ ^ E^/ujL where the formation time k of (2.9) is large compared 
to L, corresponding to situations like Fig. 7a. Recall that the LPM effect occurs when the 
relative phase K^{Xi — X2)'^ is small for space-time points Xi and X2 corresponding to 
collisions. For collisions spread out over a distance L as in Fig. 7a, this relative phase is of 
order uL{l - cos 6*) ~ luLO'^. We have LPM suppression if uLO"^ <^ 1, and the amount of 
LPM suppression of the effects of collisions with the medium turns out to be given by the 
square of this relative phase: 



I give a brief review in Appendix B 1 of why this is the amount of suppression. Putting the 
formation length (2.9) into (2.11), 



as depicted by the fall-off of e with decreasing Q± shown in Fig. 11. 

For small compared to the typical transfer of (<5x)typ ~ the solid line in Fig. 11 
deviates from (2.12), with the latter indicated by a dotted line. This qualitative difference 
will not matter to the eventual conclusions of this paper, but I will take a moment to 
explain it for the sake of completeness. In earlier discussion, I slightly oversimplified when 
asserting that bremsstrahlung is suppressed if the photon angle 6 is large compared to the net 
deflection angle A6 ~ Q_l/E of the charged particle. Consider the trajectory shown in Fig. 
12. Here, the net deflection of the trajectory is zero but the intermediate deflection is non- 
zero. The typical bremsstrahlung photon angle is then of order the intermediate deflection 
angle, which I show in more detail in Appendix B 2. Now consider cases of many multiple 
scatterings within a formation time. The rare cases where multiple scatterings produce 
smaller-than- typical deflections are dominated by situations where the first A^/2 scatterings 
produce a typical deflection, which by random chance was nearly canceled by an opposite 
deflection from the second A^/2 scatterings.^^ The dominant photon angle 6 will then be 
determined by the intermediate deflection ~ {Qi.)typ/E rather then the total deflection 
~ Q±/E that was used in (2.8) and (2.9). The upshot is that Q± <^ (Q±)typ collisions 
behave just like Q± ~ (Q±)typ collisions in terms of the medium effect on bremsstrahlung. 



For example, consider a simple one-dimensional random walk. The average displacement after N steps 
grows as N^^"^. If you look at the small subset of random walks which have zero displacement after N 
steps, their average displacement after half of those steps still grows as N^^^. 




(2.11) 



(2.12) 
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FIG. 13: Medium modification to probability for emitting a high-energy bremsstrahlung gluon of 
frequency ~ u: probabihty distribution vs. Q\. In (b), Qj. represents the transverse momentum 
picked up in a distance Zoo('^)- 



3. Putting it together 

For a thin medium, we can now find a parametric result for tlie medium effect APbrem 
on tfie probability of single-photon bremsstrahlung for photons of frequency ~ uj through 
scattering processes with total momentum transfer ~ Q^. Neglecting collincar logarithms, 
it is simply the probability of the underlying scattering event times a factor of ea for the 
associated bremsstrahlung. Multiplying (i) Fig. 1 for dPscatt/d{Q\) times (ii) Fig. 11 for e 
times (iii) a gives 

which is depicted in Fig. 13a. There is an additional logarithmic factor shown for the high- 
tail in Fig. 13a that is not included in the product (2.13). This is a collinear logarithm 
that I will explain in a moment. Recall that the notation APbrem is defined in terms of the 
frequency spectrum by (1.21). 

Fig. 13b shows the corresponding result for a thick medium L ^ /oo('^)- Since relevant 
formation lengths in this case will not exceed /oo(t^), we can break the problem up into 
independent probabilities for each section of medium of length l^. For a section of length 
/oo, we have e ~ 1. So (2.13) is modified to^^ 



C^(APbrem) L 



d{,Ql) /oo(^) 



rfP 



scatt 



d{Ql) 



X a, (2.14) 



where here Q_l refers to the transverse momentum transfer over a length of order loo{oj)- 
The parametric behaviors shown in Figs. 13a and b agree for the dividing case of L ~ /oo(<^)- 
When discussing "thick" media, 1 have implicitly assumed that a single bremsstrahlung 
analysis of the medium effect remains adequate. In particular, the media should be small 
compared to the stopping distance for the high-energy particle: L <^ L^/a. The stopping 



The dotted line in Fig. 13b showing the A''=l result for Q\ <C qloo{^) is determined by (2.13) instead of 
(2.14). In the N=l approximation, Zf can exceed loo{^) for small total deflection angle Q^jE. 
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FIG. 14: (a) initial and (b) final state collinear radiation from a single scattering with the medium. 
For comparison, (c) shows the case of colhnear radiation when there are no medium colhsions but 
only the initial hard scattering event that created the high-energy particle moving in the same final 
direction. 

distance is where AE ~ E, as can be read off from Fig. 3 for QCD or later from Fig. 17 for 
QED. 

4- Collinear logarithms 

I will now discuss collinear logarithms associated with bremsstrahlung, which were ignored 
in the previous analysis. For simplicity, consider the case u <^ E of soft bremsstrahlung, 
where the charged particle trajectory can be thought of as a classical source for the elec- 
tromagnetic field. For further simplicity, start by considering a trajectory corresponding to 
exactly one scattering from the medium, as in Fig. 8c, rather than more complicated trajec- 
tories like Figs. 8a or 7a that include multiple small-angle scatterings. In this case, there is 
a collinear logarithm in the bremsstrahlung probability associated with small photon angles 
6 <^ relative to the final particle direction, as in Fig. 14a. There is also potentially a 
similar logarithm associated with coUinearity with the initial particle direction, as in Fig. 
14b, but this second logarithm can be suppressed by the LPM effect. 

Throughout this discussion, I will assume that E and u) are large enough compared to 
effective masses that I can treat the charged particle and photon as massless. So I will not 
keep track of the cut-off of collinear logarithms due to masses. 

In the case of final-state coUinearity shown in Fig. 14a, the enhancement of the 
bremsstrahlung probability at small angles is the same as that for the vacuum process 
of Fig. 14c. There is therefore no corresponding collinear logarithm in the medium effect 
APbrem, which expresses the difference between the two. More detail is given in Appendix 
Bl. 

Now consider coUinearity with the earlier direction of the particle trajectory, as in Fig. 
14b, and let be the small angle that the photon makes with that direction. The cor- 
responding collinear logarithm will be cut off at small when the formation length (2.9) 
becomes large compared to the length of that segment of the trajectory, because then the 
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photon cannot resolve the difference between the particle trajectories of Figs. 14b and 14c. 
A generic collision in the medium will be a distance of order L from the start, and so the 
angles which contribute to a coUinear logarithm must satisfy kiOi^) < L, which is 



\^<L. (2.15) 



The angles which contribute to the coUinear logarithm are therefore 



< < A^. (2.16) 

The coUinear logarithm appears when there exists such a hierarchy of angular scales, and it 
is then 

This is the logarithmic factor shown on the large Q\ tail of Fig. 13a. The range (2.16) only 

exists if Q\ > E'^/uL. 

So far, I have considered single scattering processes like Fig. 14a-b rather than actual 
cases of interest to this paper, such as Fig. 8a. The angle that the photon makes with the 
trajectory preceding the relatively large angle collision in Fig. 8a is smeared out by multiple 
soft scatterings, which deflected the particle by an angle of order {Q±)typ/E ~ ^yqL/E^. 
The angular range (2.16) contributing to a coUinear logarithm is then replaced by 



^^\\^Ah\<0i,<A9. (2.18) 




The first case on the left-hand side dominates when L -C /oo('^), as in Fig. 13a. For L ^ 
Zoo (a;), the relevant length scale is loo{<-o) rather than L, and the logarithm becomes 

f#gl^ . ,„ . (2.19) 



(^i7)min/ 

This is argued in more detail in Appendix B3. The logarithmic factor (2.19) is the one 
shown on the large Q\ tail of Fig. 13b. 

Now return to the case of final-state radiation in a case like Fig. 8a. Once the particle 
leaves the medium, there is always a semi-infinite straight line segment of the trajectory to 
which a photon can become coUinear. So, unlike the case just considered, the effect of mul- 
tiple soft scatterings before the particle leaves the medium cannot suppress the production 
of final-state coUinear photons: photons can be produced at arbitrarily small angles (if the 
charged particle is treated as massless) by being produced after the very last scattering in 
the medium. But this contribution cancels when we subtract the vacuum contribution to 
get APbrem, just as discussed earlier for the case of Fig. 14a. 



D. Bremsstrahlung spectrum and energy loss 

1 now want to integrate over Q\ to find the spectrum APbrem = d{I — I^gj^)/dui as a 
function of uj. To see visually what values of Q±_ dominate the integration, it is useful to 
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FIG. 15: This is Fig. 13 multiplied by Q"^^. (a) assumes that l^{uj) / ^/\og <C ^ lao{uj); the figure 
is similar for L <C loo('^)/\/log, but then the right-hand (A^=l) peak would be the higher one. 



multiply the previous results of Fig. 13 by a factor of Q\ so that (i(APbrem)/'^(<5i) becomes 
the logarithmic derivative (i(APbrem)/c?(Qi)- The result is shown in Fig. 15. In the 
thin-media case, there are two peaks, corresponding to two different scales that will 
give the dominant contributions to the integral. One scale is the scale Q\ ~ qL of typical 
scattering events, corresponding to HO processes like Fig. 7a. The formation length Zf in this 
case is large compared to L. The other, larger scale corresponds to rarer, larger-angle 
scattering events that are well approximated by the A*" = 1 approximation and which have a 
shorter formation length than the typical scattering events. The scale Q\ ~ E"^ jioL of the 
right-hand peak in Fig. 15a corresponds to the case where this shorter formation length is 
of order L. The tail at yet larger corresponds to yet shorter formation lengths, as were 
depicted in Fig. 8. This tail falls off because no further gains are made in the suppression 
factor e by further decreasing /f below L, but the probability of the underlying scattering 
event decreases. 

Which of the peaks of Fig. 15a dominates in the thin-media case of L ^ ^oo(^) depends 
on the exactly how small L is. The HO and A^=l peak heights are 

APbrem(a;) - - a (-^\ ' (HO) (2.20) 



and 

respectively. The A'"=l peak dominates when 

L « ^""^""^ . (2.22) 

[ln(gL/m^)]^/' 

The (5j_-integrated result APbrem is of order the highest peak in Fig. 15. The result is 
shown vs. L in Fig. 4, provided one takes the QED formula (2.6) for Zoo = ^oo(<^) instead of 
(1.22). 

Now consider the dependence on uo when L is fixed. This behavior is shown in Fig. 16 
except that I multiply by an extra factor of uj to plot the contribution 

^^^P' " ^^brem(c^) (2.23) 
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to the medium effect on average energy loss from pliotons witli frequency ~ uj. 

Fig. 16c shows the case for thick media L ^ Loo, where L^o is given by (1.1) and 
represents the typical formation length for the case uj ^ E. In QED (unlike QCD), the 
typical formation length /oo(<-^) given by (2.6) grows with decreasing uj. For uj relatively 
large, loo{oj) will still exceed L, and APbrem will be given by the peak height 

APbrem(a;) ~ « ^ (2-24) 

of Fig. 15b. When multiplied by c<j, this gives the corresponding formula shown on the right 
of Fig. 16c. This formula works until a; gets small enough that /oo('^) ^ which occurs 
at a; ~ E"^ /qL'^. For smaller u, first the HO peak and then the A^=l peak of Fig. 15a will 
determine APbrem, giving (2.20) and (2.21) respectively, corresponding to the other formulas 
shown in Fig. 16c. 

In the case of L ^ Lqo, shown in Figs. 16a and b, we always have /oo('^) ^ and one 
or more of the stages just described are bypassed. 

Finally, integrating over a;, the total medium contribution /S.E to average energy loss just 
corresponds paramctrically to the maximum in Fig. 16. A sketch of the resulting dependence 
of AL' on medium thickness L is shown in Fig. 17. This result is qualitatively different from 
the QCD result previewed in Fig. 3, for reasons which will be explained in the next section. 
The LF' dependence of the QED result for sufficiently small L has been discussed previously 
by Peigne and Smilga [28]. 



III. BREMSSTRAHLUNG IN QCD 
A. Review of LPM Effect in QCD 

The major difference between bremsstrahlung in QCD and QED is that a bremsstrahlung 
gluon carries charge and so, like the particle that radiated, can also t-channel scatter from 
the medium. It is easier to deflect a lower momentum particle than a higher momentum 
particle and so, for the case uj <^ E, it is scattering of the gluon rather than the original 
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HO 



FIG. 17: Total medium modification to QED energy loss AE, shown vs. medium length L. 



particle that dominates determination of the angle 9 between the two. The medium effect 
on bremsstrahlung is therefore dominated by 



(3.1) 



instead of the corresponding QED angle (2.8). Correspondingly, here Q± is the transverse 
momentum that the bremsstrahlung gluon, rather than the original particle, picks up in a 
formation time. For the case to E, the deflections (3.1) and (2.8) of the gluon and the 
original particle arc parametrically the same size, up to details of group Casimirs related to 
whose Q_L we consider, which I will not bother to distinguish in my parametric estimates in 
this case. So I will use (3.1) for the entire range uj < E, assuming as always that 1—x is not 
small. 

The formation length corresponding to (2.1) and (2.9) is then 



If 



Q 



(3.2) 



The argument for the size of the formation length loo{oo) in an infinite medium then goes 
through just as in (2.4-2.6) for the QED case, but with A6' and E replaced by (A6')g and 



u, so that 



(3.3) 

V yA 

as quoted earlier in (1.22). The main qualitative difference between QCD and QED 
bremsstrahlung in a medium is that the QCD formation length /oo(^) decreases in the soft 
limit of decreasing ui due to the ease with which a soft gluon is deflected, whereas in QED 
looi'jj) increases with decreasing u. 



B. Bremsstrahlung spectrum and energy loss 

The analysis of the dependence of the bremsstrahlung problem is basically the same 
as in QED, but with the modiflcations described above concerning the formation length. 
The QCD versions of Figs. 11, 13 and 15 are given by 18-20. The only change in these 
flgures is the replacement of by a; and the clariflcation that q is qx in the case cu <^ E. 
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suppression QCD 




FIG. 18: The QCD analog of Fig. 11 for the suppression factor e. 




FIG. 19: The QCD analog of Fig. 13 for the bremsstrahlung probability dependence on Q^. 



There are group factors associated with each power of a, but I will not keep track of these 
in the figures. (Similarly, I will not distinguish between H and the density n in figures.) 
The nature of the coUinear logarithm in the QCD case is reviewed in Appendix C. 

A quick check can be made of the parametric estimate e ~ q\L^ /uo'^ shown in Fig. 18a for 
typical scatterings Q\ ~ (Qj_)typ for thin media. The medium effect on the bremsstrahlung 
probability is then of order 

Csae ^ ^ (3.4) 

by (2.7), where here I've included the factor of Cg associated with the a for the coupling of 
the bremsstrahlung gluon. Typical scatterings are described by the HO approximation, and 
(3.4) correctly reproduces the parametric dependence of the known HO result (1.8) for the 
thin media modification to the spectrum of gluon bremsstrahlung. 

Now return to the general problem. Evaluating the integration over Qj. by the peak 
heights in Fig. 20, the dependence of the medium modification APbrem on medium thickness 
L, for bremsstrahlung gluons with a frequency of order ~ a;, is given in Fig. 4, but this time 
with the QCD value (3.3) for loo{uj) instead of the QED version. 

Prom the peak heights of Fig. 20 or equivalently from the results for APbrem (<^) in Fig. 4, 
one may extract Fig. 21 showing the medium effect on the energy loss spectrum as a function 
of frequency uj. This figure is qualitatively very different from the QED version of Fig. 16 
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FIG. 20: The QCD analog of Fig. 15. For L ^ Zoo('^)/\/Tog, the right-hand peak in (a) is the 
higher one. 



d\E dAE dAE 




qL'' qLHog Eft qL^ E/2 /^.n E/2 

HO N=l HO HO 



FIG. 21: Medium modification to QCD energy loss AE: contribution of bremsstrahlung gluons 
with frequency ~ u. This is the QCD analog of Fig. 16. 



because of the qualitative difference in lao{uj). In QED, small u leads to larger formation 
lengths and so more LPM suppression, which is why the QED figure was dominated by 
Lo of order the largest scale, uj E. In QCD, small u leads to smaller formation lengths 
and so less LPM suppression, which is why for thin media the QCD figure has significant 
contributions from u <^ E. 

Figs. 21b-c show that the typical scattering processes captured by the HO approximation 
dominate (at leading-log order) not only for length L large compared to the typical infinite- 
medium formation length L^, but also in the entire range L ^ Loo/ylog, which includes 
L ~ Loo- For the case L <^ L^o/y/log of Fig. 21a, the situation is more complicated. 
Obtaining the total AE corresponds to integrating the curve in Fig. 21 with d{lnu!) = duj/uj. 
The peak in Fig. 21a gives an HO-dominatcd contribution of order the peak height, 

C.agAi-^ ~ C.CAaWL^n ( ^ ) (HO contribution). (3.5) 

The long flat plateau at larger u gives a contribution of order the height of the plateau times 
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a logarithm of its range: 



CsCAa'HL^ X — ~ CsCAa'AfL' In ( . (3.6) 

Over the frequency range gA-^^log <^ uj <^ E of this plateau, the value of u d{AE)/duj is 
well approximated by the N=l approximation, with the scattering events which dominate 
the medium effect having the form of Fig. 2b. Adding the two contributions (3.5) and (3.6) 
together gives the final result (1.18) quoted in the introduction. As described there, the HO 
contribution continues to dominate the A^=l contribution for L all the way down until the 
L.^, ~ {rri^/ q^Ey^^Loo of (1.20). This is parametrically very different from the QED case of 
Fig. 17, where HO processes dominate AE only down until L ~ L^/\/]og. 

It is important to note that there is no (A^=l)-like contribution for to <^ qxL?' in Fig. 3a 
that is separate from the HO contribution and that would parametrically reduce the lower 
scale qAL"^ in the logarithm of (3.6) to make the sum (1.18) of (3.5) and (3.6) larger. To 
check, 1 will estimate the contribution 5{/S.E) to /\E from (A^=l)-like scattering (Fig. 2b) 
with uj < r/A^v^, as opposed to HO- like scattering (Fig. 2a). The range uj <C ^a-^^ in Fig- 21 
corresponds to L looij^) and so to Fig. 20b. The contribution to APbrcm from the A^=l 
part of that curve is of order the height in Fig. 20b where the HO and A'^^l curves meet: 

\n[qAL/mf^) l^{uj) \n{qAL / mf^) V oj 

Multiplying by u and integrating over the range uj < (IaL^ under discussion gives an addi- 
tional N=l contribution to AE of 

5{AE) L 5(APbre..)] , ~ j7§T7^ ^ C.CacW. (3.8) 

This contribution from (iV=l)-like events with u) < qA^^ is sub-leading in logarithms com- 
pared to the HO contribution (3.5) and the total result (1.18), and so it can be ignored in 
a leading-log analysis. 



IV. COMPARISON TO ZAKHAROV'S ANALYSIS 
A. The Puzzle 

In Ref. [20], Zakharov argued that the HO approximation should be expected to break 
down when L is less than or order the infinite-volume formation time. In this paper, I 
have argued that the HO approximation does a little better than that if one consistently 
treats logarithms as large. For the case of the medium modification to the bremsstrahlung 
spectrum, depicted in Fig. 20, the HO approximation dominates as long as L ^ /oo(^) / Vlog, 
which includes L ~ /oo(^)- In this section, I will paraphrase Zakharov's argument and resolve 
the slight difference in conclusion. 

First, I need to briefiy review the formahsm for doing a full calculation of the gluon 
bremsstrahlung spectrum [3], which was originally developed for finite media by Baier, 
Dokshitzer, Mueller, Peigne, and Schiff [4-6, 14] and by Zakharov [7, 8] . A brief summary in 
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my own notation, which is suited for discussing problems where the particles in the medium 
are not fixed scatterers, can be found in Ref. [9]. The spectrum is given by [9] 



oo 



dw^ [x{l-x)EY Jo ^ 

VB,-'VB,{G{B2,t2;Bi,ti)-G,^,{B2,t2]Bi,ti)}] , (4.1) 

where G{B2,t2; Bi,ti) is the Green's function for a two-dimensional quantum mechanics 
problem with the time-dependent, non-hermitian Hamiltonian^^ 

H{t)=6E{pB)-tT,{B,t). (4.2) 

Here 6E describes the energy difference {Es^p+Eg^iS)—Es,p+k between (i) a high-energy parton 
of momentum P = p+k and energy E = P and (ii) the same parton with momentum p plus 
a bremsstrahlung gluon with momentum k. In the high-energy limit, if we ignore masses, it 
can be written as 

where Pb = {pk± — kp±)/P is the transverse momentum conjugate to the B of (4.2) and 
M = x{l — x)P is the "mass" of the two-dimensional Schrodinger problem: 

H{t)..^-iT,{B,t). (4.4) 

For fixed x, this Pb is proportional to the angle between k and p. The second term in (4.4) 
is 

r3(S, t) = iC-A f2{B, t) + (C, - iCA) HxB, t) + iC-A r2((l - x)B, t) , (4.5) 

where 

f,(6,t)^ |dV^(l-e^'"--). (4.6) 

Here f ei is defined by Fei = C^f ei. That is, it is the elastic scattering rate without the group 
factor Cr associated with the particle being scattered. 

The high-energy limit corresponds to large mass M in the two-dimensional Hamiltonian 
(4.4) and so will be determined by the small B behavior of the "potential" —iT^^Bji). 
Naively, (4.6) for small b gives 

f2(6,^)^|6^ (4.7) 

where Crq is formally the average momentum transfer per unit length rather than the typical 
transfer used throughout this paper. The actual small b behavior of r2 is proportional to 
6^ ln(m^6^), not 6^, which is reflected by the UV divergence of the integral (1.11) for average 
transverse momentum transfer. Cutting off this divergence by replacing Crq by the typical 
momentum transfer per unit length, as in (1.13b), corresponds to the harmonic oscillator 
approximation, so named because of the form of (4.7). 



Zakharov uses the letter p for what I call B. For a complete translation table, see the appendix of Ref. 
[9]. 
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In contrast, another analytic approach to solving the problem is to keep the full origi- 
nal form of (4.4) and instead do perturbation theory in powers of the (imaginary-valued) 
potential —iT^. This is the formal version of the opacity expansion. 

Alternatively, both approximations can be made. Consider the case of the brick problem. 
If one first makes the HO approximation (4.7) and then makes the opacity expansion, the 
opacity expansion is simply the expansion of the HO result (1.5) in powers of |a;Q| oc g oc 



-ir. 



d 



4 r4 



-^vac)HO = -xPs^g{x) [^\uJo\''L 

Parametrically, this expansion has the form 



17 I, , 187-8 I ] 



(4.8) 



^Iw^^ ~ -^vac)HO ~ 



Cm 



C,a 



Cm 



# 



# 



# 



H-#(-^^ 



UJ 



+ 



+ # 



CAQiWLMog 



+ 



+# 



CAQiWLMog 



+ 



(4.9) 



with the logarithms defined as in (1.16). The condition for the perturbative expansion of 
(4.7) to be useful is that successive terms get smaller and smaller, which parametrically is 
the condition that L ^ /oo(^)- As noted by Zakharov [20], there is no first-order term (no 
term proportional to and so proportional to the interaction Fs) in the expansion (4.8). 
But if L -C looij^) so that a perturbative treatment of the quantum mechanical problem is 
valid, then why not forgo the HO approximation and just use the full, original potential 
—iV^,. At first order, one then obtains the A'"=l result (1.9), so that 



Cm 



if \- 



(4.10) 



The fact that a perturbative expansion of the HO result should work whenever L <^ /oo(<-^), 
yet the HO approximation is clearly missing the first-order term in this limit, makes it seem 
like the HO approximation must be untrustworthy whenever L <^ ioo(<^)- 



B. Reconciliation 

The absence of the first-order term in the expansion of the HO result can be illuminated 
if one separates out from Eq. (4.1) the step of taking the real part. The origin of the HO 
result (1.5) is actually 

d ot 
UJ —{I - /vac)HO = -xPs-.^{x) Rc [In cos(a;o-L)] . (4-11) 
au) TT 

Correspondingly, using lJq = —i\ujo\'^, the perturbative expansion is 

OJ -^{I - /va.)HO = ^xPs^^ix) Re [||a;o|'L2 + ^|a;o|^L^ + ^IcuqI'L^ - ^\ooofL^ + ■■■], 

(4.12) 
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which shows the missing odd terms in the expansion. This corresponds to 



LO -^(/-/vac)HO ~ Cga Re 



+ 



(4.13) 



In contrast, the full A'"=l perturbative calculation turns out to give 



uj^{I - /vac)Ar=i - L"^ X Ps^gix) Re 



d 

djjj ' 



to 



•1 / 



L X Ps^g{x) Re 



i In 



+ 



TT 



(4.14) 



in the limit of large logarithms. Taking the real part and the small x limit, this reproduces 
(1.9). Parametrically, 



uj^{I - Ivac)N=i ~ CgaRe 
duj 



~ CsOi Re 



In 



^4 



to 



(4.16) 



Comparing (4.13) and (4.15) before taking the real part, note that the first-order terms 
are parametrically the same except that the arguments of the large logarithms are different. 
When the real part is taken, however, nothing survives of the first term in the expansion 
(4.13) of the HO result, but a term sub-leading in large logarithms survives from the A^=l 
result. The moral is that the structure of the A^=l and HO results for thin media are not 
very different before one takes the real part. After the real part, a calculation which included 
contributions from both HO and A'"=l physics, as described in this paper, would be expected 
to produce a result of the form 



'^-^i.I - -^vac)HO ~ CgO. 

du! 



(4.16) 



This is consistent with the L < loo{'^) behavior of APbrem found in this paper and shown in 
Fig. 4. 



C. Some Differences 

It is important to note, however, that a perturbative calculation to second order in —iV^ 
would not give precisely (4.16) with the HO logarithm (1.16), 

When discussing the perturbative expansion (4.9) of the HO result, I first made the HO 
approximation (4.7) and treated g as a constant given by (1.13b). Only then did I expand in 
powers of —iV^. If I instead forgo the HO approximation, then the expansion in —iT^ (the 
opacity expansion) is equivalent to an expansion in powers of the medium density ~ Hi if 
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for this purpose I treat the Debye screening mass mo as a variable independent from A/". 
However, (4.16) is not a simple power series expansion in A/", because there is a factor of 
M inside the argument to the logarithm (4.17). The second-order HO term in (4.16) can 
therefore only arise from a resummation of many terms of the opacity expansion. 

In this paper, I will not attempt to explore in detail how physics associated with the HO 
approximation can be seen to arise from resummation of terms in the opacity expansion. 
But I hope that the discussion of this section gives some insight into how the earlier results 
of this paper, based on more physical arguments, can be consistent formally with the small 
L expansions of the A^=l and HO approximations. 
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APPENDIX A: THE FORM OF dPscatt/d(Qi) 

The qualitative form of Fig. 1 — a diffusion peak from typical scatterings with a single- 
scattering tail due to rare scatterings — has a very long history. For example, a simplified 
version of MoUiere's theory of multiple scattering was given by Bethe in 1953 [30].^^ In the 
current context of leading-log approximations to jet Qj, broadening in quark-gluon plasmas, 
it has been addressed previously in Sec. 3.1 of Ref. [6] and is nicely reviewed in Appendix A 
of Ref. [28] for the particular model of scattering where dTc\/d{q\) is taken to be proportional 
to (g^ -|- m|,)~^. For the sake of completeness, I will review here the important elements 
for the current work in a model-independent way. Specifically, the important points for my 
argument in this paper are that (i) dPscatt / d{Q\) ~ for Q\ <S qL (the peak height in 
Fig. 1), and (ii) dPsca.tt / d{Q\) ~ a'^nL/Q]_ for Q\ 3> qL (the form of the large Q± tail). 

The transition between these two behaviors occurring for qL < Q\ < qL In(log) in Fig. 1 
is interesting but unimportant to my conclusions. In this paper, I treat logarithms as large 
but I treat logarithms of logarithms In (log) as 0(1). So Q\ :$> qL refers to the tail of Fig. 
1 and not to any part of the transition region. 

Though not needed for the present work, I will also provide reference to a rigorous 
mathematical generalization of the central limit theorem which demonstrates that the soft- 
scattering peak indeed approaches a Gaussian form in the limit of a large number of collisions. 



Readers may wonder how a logarithm of the form \n{cN') could possibly arise from any power series in 
Af, since ln{chf) is not expandable as a power series. Keep in mind that the form (4.17) of the logarithm 
is only meant to be valid in the limit that the argument of the logarithm is large. So, as an example, 
ln(l -|- cJ\f) ~ \n{cJ\f) when the argument is large, but ln(l + cAf) has a series expansion in Af. 
For other references, see Sec. 27.3 of the 2008 Review of Particle Physics [31]. 
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1. Review of general multiple scattering formula 



Let 

r ( ,\ C^i^scatt 1 ^-Pscatt / A i \ 

d?Q±_ TT d{Q\) 

be the two-dimensional probability distribution of Q_l at time t. I follow the standard 
development of multiple scattering by writing the evolution equation for /, which is 

dtf{Q±,t)=[ p{q±)[f{Q±-q±,t)-f{Q^,t)] (A2) 
in a uniform medium, where 

is the two-dimensional probability density for acquiring a transverse momentum kick of q± 
in a single, individual collision. The first term on the right-hand side of (A2) is a gain 
term, representing momentum change from Q±—q± to Q±. The second term is a loss term, 
representing change from Q± to Q±+q±. Take the initial condition f{Q±,0) = 5^^-*(Q±). 
The equation is solved by Fourier transformating to 

dJ{b,t) = [~p{b)-Pmfib,t) (A4) 

with initial condition f{b, 0) = 1 and solution 

/(6,t) = e-[^W-^('')l*. (A5) 

Fourier transforming back, 

p(0) - p{b) = ^^1^(1 - e^'-'^). (A7) 

2. Examples 



where 



As an example, consider the weak-coupling result [32]^^ 

dr^ CRg^Tml 

^ {2T:Yql{,ql + ml) 



(A8) 



which holds for q^^ <^ T. The result is not significantly different for q± ^ T [24, 25], so 
I shall take it as an example for discussing the entire range of q±. The behavior for 



For a brief overview in the notation used here, see Sec. II A of Ref. [25] . 
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q± -C mo is due to magnetic scattering, which is not completely screened by the Debye 
effect. The Fourier transform of (A8) gives [33] 



m - m 



Crq^T 
27r 



Ko{mBb) - In 



Note that for -C 1 this becomes 



p(0)-p(6)^%^(mD6f 



Stt 



In 



mj)b 



7e + 1 



(A9) 



(AlO) 



and for ^ 1 it is 



p(0) - p{b) ~ 



27r 



In 



(— 

\mT)b 



+ 7e 



+00 as 6 — > oo. 



(All) 



Alternatively, consider a popular model used in this subject, which is 

dr^ Crq^N 



(A12) 



where H is the number density weighted by appropriate group factors.^^ Fourier transform- 
ing, one finds 

p{b) = bKi{mBb), 



so that 



For mnb ^ 1, 



47rmr 



p(0) - p{b) = [1 - mD6i^i(mD6)] 

47rm 



D 



m - m 



Stt 



In 



(. 



and for ^ 1 it is 



p(0)-p(6)~p(0) 



7e + 



as 6 — > oo. 



(A13) 
(A14) 

(A15) 
(A16) 



The 6^ ln(mD&) at small b in both examples is a universal result of having a drei/(fq± that 
falls as l/o'l at large q±. This in turn is a universal feature of point-particle scattering that 
is Coulomb at short distances. Note that the two formulas (AlO) and (A15) have the same 
size paramctrically, so it does not matter which we use if we are interested in a parametric 
analysis. In general. 



p{0)-p{b)^CngW In 



m-ob 



for mjib <^ 1. 



(A17) 



In detail, I am using the notation of Ref. [25]. 

Eq. A15 is the correct formula at high enough energy that typical individual scatterings have q± ^ T. 
See, for example, the discussion in Ref. [25] . 
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In the second example, the finite b^oo hmit (A16) means that the final Fourier transform 
(A6) used to obtain dPsca.tt/(PQi_ contains a 5-function singularity. This can be isolated by 
rewriting (A6) as 

dPscatt _ -p(0)t 



= e + ^ [e-\m-~Pm _ e-p{o)t] g-^^>•Qx_ (a18) 



The coefficient of the 5-function is just the probability exp(— Fcit) that there are no scatter- 
ings whatsoever. For large Fdt this can be ignored. This 5-function does not appear in the 
first example (A8) since in that case Fd = /^^ d^Vd/ d^q^ = oo because of the l/q1_ infrared 
behavior of (AS) due to magnetic scattering.^^ Either way, I will ignore the 6 function in 
the rest of this discussion. 



3. The height of the peak 

For the height of the peak if Fig. 1, we just need to evaluate the regular (i.e. non-5- 
function) term in (A18) at Q±=0: 

For large enough t, this integral is dominated by small b determined by [p(0) — p{b)]t ~ 1. 
Using (A17), this is 

CRg^Ntb'^lTL f ^ ) ~ 1, (A20) 



and so 

1,2 



\mY)b 

ml 



(A21) 

qRt ^ ' 



The condition rriDb <^ 1 that I have used is then satisfied provided t ^ ml/ {Crq^N), which, 
neglecting group factors, is of order Ijg^T. This is just the condition that there are many 
g_L~mD scatterings, which I have assumed throughout this paper. Prom (A21), the size of 
the integral (A19) is then 

Replacing t by L, this is just the peak height depicted in Fig. 1. 

The Q±—0 result will be a good approximation whenever the exp(i6- factor in (A18) 
is approximately 1 up to and including b values of order (A21). So dPscatt/ d'^Qi, ~ I/qkL 
for Q± <^ b~^ ~ yJq^L^ just as shown in Fig. 1. For larger Qx, the oscillating factor will 
cause dPscatt/d'^Q± to fall. For much larger Q±, the oscillating factor exp{iQ± ■ b±) causes 
the integral in (A18) to be dominated by even smaller b, in which case one may approximate 

g-[p(o)-p(6)]i ^ _ j~^Q^ _ ^^23) 



This infrared divergence will be cut off at q± ^ g^T by non-perturbative physics. The resulting (5-function 
term will in any case be exponentially small when t S> I/ttid, which has been assumed throughout this 
paper. 
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Completing the Fourier transform, this just gives that the large Q± behavior of dr/(PQ± is 
given by the single-scattering formula dVei/d^qi, with = Q±, producing the tail of Fig. 
1. 



4. The Gaussian shape of the peak 

Ref. [6] discusses the Gaussian shape of the peak by noting that ln(6^) is a slowly varying 
function and so making the approximation of replacing this logarithm by a constant in (A18). 
This is the harmonic oscillator approximation and gives a Gaussian result for rf-Pscatt/c^^Q±- 
Ref. [6] notes that this approximation breaks down for the tail, which is generated by the 
non-analyticity of In 6 at 6=0. Some readers may wonder, however, if the argument that the 
shape approaches Gaussian can be made more rigorous. If the distribution d^Vd/ cpqi_ of 
single-scattering rates had a finite variance (g^) < cxo, then the the approach to a Gaussian 
shape would be guaranteed by the central limit theorem. Finite variance is a sufficient 
condition for the central limit theorem, but it is not a necessary condition. A necessary 
and sufficient condition may be found in Theorem 8.1.3 of Ref. [34].^^ I will simplify their 
condition to the radially-symmetric case of interest hcrc.^^ Using the notation of this paper, 
the condition can be written in terms of the truncated variance 

gA= / (qx-n)V(g±)^(A-|q±-n|)), (A24) 

where 6{z) is the step function, could be a vector of any dimension, p{q±) is a probability 
density in that vector space, A is a cut-off, and n is a unit vector in any direction. The 



The limit is non-uniform, which means that for finite t there wiU stiU be non-Gaussian tails. However, as 
t increases, the region of over which Gaussian is a good approximation becomes larger and larger in 
units of the width of that Gaussian. 

Ref. [34] gives a condition for distributions of vectors. The specialization to one dimensional distributions 
has an older history: see Theorem 1 on p. 172 of Ref. [35] and references therein. Also, the theorem is 
formulated for a sum of a finite number of vectors drawn from the distribution p. In our case, one should 
think of p as the probability distribution for picking up transverse momentum q± in a very small time 
interval At, and the total Q±_ is then the sum of these transfers, with the limit ^ taken at the end 
of the day. 

Ref. [34] , which also applies in the absence of radial symmetry, defines F{x) = J^^ (x ■ <7_l)^ p{q±) ^'(1 ~ \x ■ 
q±\). Their condition is that there exists (i) a function f(t) from to the set of linear transformations 
on R'^ with f{Xt) f{t)^^ X^^I (where / is the identity operator) as t ^ oo for some E > and for all 
A > and (ii) another function R{t) from R+ R+ with R{\t)/R{t) A^^, such that (iii) 

lim^ML^=0(.) 
t^oo R{t) ' 

for some (p{x) > whenever ^ a; in M** - {0}. See pp. 96, 125-7, and 293 of Ref. [34]. For the 
radially symmetric case, F{x) is proportional to |a;p times what I have called ^a, and their x is my 
n/A. Workable choices for the specific isotropic case of interest here, where oc InA at large A and so 
F{x) oc ln(l/x) at small x, are /(t) = tl and R(t) = Int. 
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necessary and sufficient condition for tfie central limit theorem is that 

lim ^ = 1 (A25) 
A^oo q\ 

for all A > 0. This condition applies to the case relevant in this paper, where oc InA in 
the limit of large A. 



APPENDIX B: MORE DETAILS ON SUPPRESSION e IN QED 

1. Review: Why e ^ {K • Ax)^ in QED 

Here I will briefly review the parametric formula (2.11) for the suppression factor e. For 
simplicity, I will focus on QED in the soft bremsstrahlung limit, u <^ E. In this limit, 
the charged particle trajectory can be thought of as a fixed, classical source J'^{t, x) for 
the electromagnetic field, and the bremsstrahlung amplitude is proportional to the Fourier 
transform J^{u!, k). For further simplicity, I will focus on a comparison of a single scattering 
from the medium, as in Figs. 14a-b, with the vacuum case of Fig. 14c. 

The Fourier transform J^[K) of the current for the trajectory shown in Figs. 14a-b is 



K-V{ 



JK-Xi 



(Bl) 



where K = {uj,k) = {k,k) is the photon 4-momcntum, and V = {l,v) with the charged 
particle velocity v before (Vi) and after (Vf) the medium collision. The medium effect on 
the bremsstrahlung probability is then 



\e-Jr-\e 



7r oc 

Ivac ^ 



s-V, 



-iKAX 



-1) + 





to 


e-Vi 


K -Vf 




K-V{ 



(B2) 



where e'^ is the photon polarization. The first term in (B2) corresponds to Figs. 14a-b and 
the second to Fig. 14c. Expanding the first square and summing over polarizations gives 



J2 [k • J\' - k • J\L] oc 2 [1 - cos(X • AX)] J2 Re 



e*-Vi f e-Vi e-Vi 



K-Vi\K-Vi K-Vi 



(B3) 



For K ■ AX <^ 1, the bremsstrahlung probability is therefore proportional to {K ■ AXy, as 
in (2.11). 

The e ~ (AT ■ AX)^ suppression factor favors larger photon angles relative to Aa; over 
smaller ones. So it's important to review what sets the upper limit on photon angle in the 
current context. The difference (B3) between the bremsstrahlung probabilities in medium 
and in vacuum will sometimes be positive and sometimes negative. Consider the average 
of (B3) over rotations of the direction of Vf around the axis defined by Vi. Let and 9fi 
be the angles that k and V{ makes with Vi, and let 0fi represent the azimuthal angle being 
averaged over. Ignore mass effects, so that Vf is a unit vector. The averaging gives 



e-Vf 
K-Vf 



4>ti 



e-Vi 
K-Vi 



if 9^1 > 9ii, 



(B4) 
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V, 



FIG. 22: An isolated collision. 



where (■ ■ ■ )<^,fi indicates averaging over 0fi. So the corresponding average of the medium 
effect (B3) vanishes if is larger than the deffection angle A9 = 9{i of the charged particle 
trajectory in Fig. 14. For multiple scatterings involving a number of consecutive trajectory 
directions Vi, V2, Vtv = v^, one can start by averaging over rotations of Vj^ about fN-i 
and then work back recursively to show that a similar cancellation occurs if 9^^n-i > ^n,n-i 
for all n. 

Eqs. (B2) and (B3) also show the behavior associated with coUinear logs discussed in 
Section II C 4. Compare (B3) to the usual calculation of bremsstrahlung from an isolated 
collision like Fig. 22: 

e-Vi e-Vf " 



JP oc 



K-v. ■ K-v, - <^«' 

The coUinear divergences associated with photons coUinear with the initial or final direction 
correspond to the divergence of (B5) as K -Vi or K -Vf vanishes, respectively. If polarizations 
are summed over, this corresponds to divergences like 1/9?^ and i/9j^, where 9i^ and 9iy 
are the corresponding angles. Integration over angles (PD, then gives the usual coUinear 
logarithmic divergences. In (B2), however, the cancellation with the vacuum term softens 
the K ■ Vf ^ behavior, so that there is no corresponding logarithmic divergence in the 
angular integration. The other divergence, K-Vi ^ 0, is cut off once K-AX oc K-Vi becomes 
smaU enough that 1 - cos{K ■ AX) ~ {K ■ AX)'^ in (B3): that is, when K ■ AX < 1. 

2. The case Q± < (Qi)typ 

The diagram of Fig. 12 corresponds to 

^ K-v, ^ K-V^ K-V, 

^ ^iK.X, f Jf V/^\ ^ -iK.AX _ (gg) 

\K-Vi K-Vj^ ^ ^ 

where Vi is the intermediate particle direction and the initial and final directions are equal: 
V = V{. So 



JPoc 



s-V e-Vi 
+ 



le-'^-^'' -if . (B7) 



K-Vi K-Vi 

The factor \e-'^-^^ - Ip gives e - • AX) LPM suppression when K - AX < 1. The 
other factor in (B7) looks just hke the result (B5) that the individual coUisions would each 
give if they were isolated. The two terms in this factor approximately cancel each other only 
when the photon angle is large compared to the angle between V\ = V{ and Vi. In order of 
magnitude, the characteristic angle 9 of bremsstrahlung from two collisions with canceling 
defiections is therefore similar to that of two collisions with deflections in the same direction. 
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3. CoUinear logarithms for thick media 



Here I will flesh out the argument for the coUinear logarithm (2.19) that appears for 
the high-(5± tail in Fig. 13b in the thick medium case L » ^oo(<^)- Return to Eq. (2.18). 
Multiple scatterings preceding the rare, larger-than-typical angle scattering will only affect 
the bremsstrahlung process if they occur within the corresponding mean free time li{9i^) ~ 
l/cuOf^. So, the condition (2.19) should more accurately be written 



For L ^ the lower bound on which dominates is 



''^^'^^^^J^<0,,. (B9) 



This is equivalent to 



1/4 

< e,,. (BIO) 



uE'^ 

Using (2.6), the constraint (B8) can in this case be written in the form 

i^^^^^<^i.<A^^% for L»U'^). (Bll) 
A significant range exists when ql^{uj), and the corresponding logarithm is (2.19). 

APPENDIX C: COLLINEAR LOGARITHMS IN QCD 

In the case of QCD, the gluon scatters from the medium. If we again focus on the case 
of a single significant scattering, then Fig. 23 needs to be added to the situation considered 
for QED in Figs. 14a-b. For simplicity, neglect the original creation of the particle at the 
left-hand side of these diagrams and instead considers it to come from infinity. This is the 
gluon bremsstrahlung situation analyzed long ago by Gunion and Bertsch [36]. Schemati- 
cally, the amplitude for bremsstrahlung compared to the amplitude for scattering without 
bremsstrahlung is (in the high energy limit) proportional to ^'^'^^^^ 



{k^-xqA_)-£^{l-x) -£±(1 -x) ^fei (fc±-g±)-£±(l-a;) 

2 



|fe_L — xq_Lp k, |fc_L — q_Lp 



(CI) 

where _L is defined relative to the initial particle direction, and a characterize the trans- 
verse momentum transfer and adjoint color index associated with the collision, and fc^ and 
h characterize the transverse momentum and adjoint color index of the final bremsstrahlung 
gluon. In the corresponding QED calculation, the first two terms cancel in the limit that the 
photon angle is large compared to the deflection angle of the charged particle {k±^ ^ xqi_). 
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My qj_ and fej, are respectively the l±_ and q±_ of Ref. [36] . 



37 



FIG. 23: A QCD addition to the diagrams of Fig. 14a-b. 



In QCD, they do not because the color generators do not commute. Instead, in the soft gluon 
case, the cancellation occurs between all three terms of (CI), when the bremsstrahlung gluon 
angle becomes large compared to the deflection of the gluon due to q± {k± ^ q^). That is, 
in the notation of (3.1), bremsstrahlung gluons associated with the collision decouple when 
e > (A^)g. 

The amplitude proportional to (CI) diverges when any of the denominators goes to zero, 
leading to collinear divergences. Now return from the Gunion and Bertsch problem back 
to the case of Figs. 14 and 23 where the original particle is created at the beginning of the 
trajectories shown. The collinear divergence corresponding to Fig. 14a, i.e. the square of the 
first term in (CI), will cancel against the corresponding virtual correction to bremsstrahlung 
in the vacuum case. The rest will give a collinear logarithm that will be cut off at small 
angles when the LPM formation time becomes of order L, just as in the QED case. One 
way to understand this is to note that collinear divergences correspond to the intermediate 
particle states in Figs. 14a-b and 23 going on shell. However, for the case of Figs. 14b and 23, 
the intermediate particle lines can have length at most L, which introduces an uncertainty 
in their energy of order 1/L. This provides a lower bound to how on-shell those particles 
can be and so cuts off the corresponding collinear divergences. 



APPENDIX D: CONTRIBUTION TO A^; FROM x^l 

The HO formula (1.2) for AE comes from taking the small x approximation to the 
spectrum (1.5): 



d 2Csa 

^ - JvacjHO In 

duj TX 



cos ■ 



2uj 

In the thin- media limit, integration over uj is dominated by small x of order 



(Dl) 



(D2) 

Integrating (Dl) over uj from zero to infinity gives (1.2). But this procedure implicitly ignores 
the possibility of an additional contribution from small 1—x in the case of bremsstrahlung 
from a quark (q— ^-gq) or anti-quark. 
Define 

d 

AE^ = j duuo—{I-U,,) (D3) 
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to be the average energy loss of the quark in q^gq. In the thin media hmit, one contribution 
to AEq is the small x approximation just discussed. But there is another contribution from 
small 1—x. In the limit of small 1—x, the spectrum (1.5) becomes 



dcu TT 



In 



cos . 



—iqpL"^ 



2{l-x)E 



(D4) 



The u) integral of (D4) is dominated by small 1—x of order 



1 — x 



E ' 



(D5) 



and the integration gives an small 1—x contribution of ^C^aq-pL'^ to A£^q. Adding this to 
the small-x contribution of (1.2), 



(A£;)q,HO - \CMqK + \qF)L'. 



(D6) 



Now instead consider the average energy loss of the leading parton for q^gq, which is 

■.E/2 ^ rE 

duJU—{I - /vac) + 

'E/2 



AE^iead = / dLLlLJ—{I- /yac) + / du {E - Lj) —{I - /vac)- 

Jo duj Je/2 duj 



(D7) 



The replacement of the uj factor by E—uj in the second term produces an additional suppres- 
sion of (D5) to the small 1—x contribution. As a result, the small x contribution dominates 
in the thin media limit, and so A/?icad is simply the /S.E quoted in the main text. 

If one is interested in understanding the parametric dependence of the q— >gq 
bremsstrahlung spectrum for the case of small 1— a;, it is easy to adapt the QCD results 
of the main text. In this case, the final-state quark is the particle that is most easily scat- 
tered and so the one whose scattering sets the scale of the LPM effect. Paramctrically, the 
result for the bremsstrahlung probability when u > E/2 will look just like Fig. 4 but with 
the replacement 



'E-uj 

qp 



(D8) 



instead of (1.22). 
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